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. ^ . 1 INTRODUCTION 



ABSTRACT 

We investigate the global dynamic stability of spherical clouds in the Broad Line Region 
(BLR) of Active Galactic Nuclei (AGN), exposed to radial radiation pressure, gravity of 
the central black hole (BH), and centrifugal forces assuming the clouds adapt their size ac- 
cording to the local pressure. We consider both, isotropic and anisotropic light sources. In 
both cases, stable orbits exist also for very sub-Keplerian rotation for which the radiation 
pressure contributes substantially to the force budget. We demonstrate that highly eccentric, 
very sub-Keplerian s table orbits may be found. This gives further support for the model of 
iMarconi et al.l (l2008l) . who pointed out that black hole masses might be significantly underes- 
timated if radiation pressure is neglected. That model improved the agreement between black 
hole masses derived in certain active galaxies based on BLR dynamics, and black hole masses 
derived by other means in other galaxies by inclusion of a luminosity dependent term. For 
anisotropic illumination, energy is conserved for averages over long time intervals, only, but 
not for individual orbits. This leads to Rosetta orbits that are systematically less extended in 
the direction of maximum radiation force. Initially isotropic relatively low column density 
systems would therefore turn into a disk when an anisotropic AGN is switched on. 
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Recently, the BLR has received much attention, because the line 
widths, in combination with the size of the region, measured for 
example by reverberation mapping, allows a de termination of the 
mass of the c entral super-massive black hole (e.g. lBentz et alj2009l ; 
lNetzeiil2009n . which is essential if one wants to place AGN in the 
context of general galaxy evolution. From such studies, one finds 
that AGN generally have normal sized black holes, as expected for 
the size of their host galaxy, and that most galaxies must have been 
active for one or more times in their life. The measurement of the 
black hole mass is affected by the relative importance of radiation 
pressure and gravity, which should dominate the force budget. This 
is connected to the properties of indi vidual BLR clouds a nd the 
dynamical configuration of the system JMarconi et al.ll2008|). 

The BLR i s a standard A GN ingredient (e.g. lOsterbrockl 19881 ; 
IPetersonI [19971 ; lNetzeill2008l) ; It is located inside the obscuring 
torus, above and below the accretion disk, and at a distance of 
fractions of a parsec from the super-massive black hole. The line 
emission is powered by photoionisation by a broad-band contin- 
uum due to the innermost parts of the accretion disk. Photoion- 
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isation also governs the thermodynamics of these clouds; They 
have a stable equilibrium temperature of order 10 K, indepen- 
dent of their locations. The critical densities of observed and sup- 
pressed emission lines and pressure equilibrium considerations 
JFerland & Elitzuii Il9841) constrain the number densities in the 
clouds to Hci = 10 cm^ . Signatures from partially ionised 

zones are typically observed. Therefore, the Stromgren depth is 
a lower limit of the cloud size. An upper limit may be found by 
requiring the lines to be smooth, which demands a certain min- 
imum on the number of individ ual entities. This results in cloud 
sizes of about r^ = 10 cm JLaor et afluOOg . and references 

therein). In agreement with this, detailed photoionisation models 
I IKwan & Kr olik 1981) constrain the column densities to A'd > 
10 cm^ . Because BLR clouds are optically thick, high ioni- 
sation lines are produced only by the illuminated surface of the 
clouds. In such lines, one observes generally the far side of the 
BLR, wherefore outflows manifest themselves as redshifts, inflows 
as blueshifts. 

Much information has also been gathered on the dynamical 
state of the BLR. Inclination matters; the statistics of line widths ex- 
cludes Keplerian rotation in a flat disk, but is instead consistent with 
a thick disk configuration with v /g ~ 2.5, where v is the rotational 
velocity and CT the turbulent one JOsterbrockll 19781) . For radio-loud 
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AGN, the ratio of core to lobe power correlates with the width of the 
broad lines in the sense that the systems ob served at a line-of- sigh t 
close to the jet axis have the na rrower lines (Wills & Bro wnell98a) . 
This finding was confirmed bv lJarvis & McLure ( 2006 ). They find 
that the average line width of flat spectrum radio sources, which 
are seen more face on, have narrower emission lines by about 30% 
than steep spectrum radio sources. Spectropolarimetry of type 2 
objects, where the BLR is hidden by a dusty torus, has revealed 
hidden BLRs in pol arised light, which established the orientation 
unification scenario dAntonuccill 19931) . In these objects, the BLR 
emission is scattered by a polar region above the BLR. The polari- 
sati on angle is perpend icular to the system axis, defined by the ra- 
dio l lSmith et alj|2004l) or UV/blue dBorguet et allbOOSh emission. 
The same authors find a preferentially parallel orientation for type 1 
objects, which provide an unobscured view towards the BLR. This 
finding is explained by an equatorial scattering region, much closer 
to the BLR. Where observed (type 1 objects), the latter dominates 
the polarised emission. The equatorial scatterer is in fact very close 
to the BLR and provid es unambiguous ev idence for a disk-like con- 
figuration of the BLR dSmith et al.ll2005l) : Approaching and reced- 
ing parts of the BLR-disk are seen at slightly different angles by the 
equatorial scatterer. This produces a noticeable difference in polar- 
isation angle for the blue and red par ts of the line . This , and o ther 
characteristics, has been observed bv lSmith et alj 1 120041 120051) for 
a sample of Seyfert galaxies. They find a continuum of polarisa- 
tion properties: Zero polarisation for near pole on objects, where 
the different contributions by equatorial scatterers in various di- 
rections cancel each other; wavelength dependent parallel polari- 
sation for intermediate inclination type 1 objects, where parts of 
the BLR and the equatorial scattering region are attenuated by the 
obscuring torus; and perpendicular polarisation for type 2 objects. 
Other evidence pointing to a disk-like nature of the BLR comes 
from the existence of objects with double peaked broad lines in- 
dicative of Keplerian motion dEracleous & Halpemll2003l 12% in 
their sample), while the majority of single peaked objects allows 
for a hidden disk component, if the inclination is small, and th e 
contribution of the flattened part is not too strong dBon et al.l2009f) . 
These findings suggest that the dominant contribution to the BLR 
kinematics is Keplerian rotation, followed by a turbulent compo- 
nent. There is also evidence for radial motion: There are examples 
in t he literature for bulk outflows, as measured by spectropolarime- 
try dVoung et al.ll2007l . and references therein). Velocity resolved 
reverberation mapping confirms the dominant Keplerian motions, 
but additionally finds evidence for bulk inflow and outflow in in- 
dividual objects (e . g. Done & Krolik 1996; Ulrich & Horn^ | 19961 ; 
Kollatschnv 20031; [Bentz et all l2008l |2009| ; IPennev et alTfeoogh . 



GaskelL ^2009.) interprets the combined evidence of velocity re- 



solved reverberation mapping as strong evidence that the radial part 
of the kinematics is dominated by inflow. 

Radiation pressure due to the central parts of the accre- 
tion disk is generally assumed to be significant in the BLR (e.g . 



Blumenthal & Mathewsll"975l ; lMathewsll986l ; lMarconi et al.l2008l ; 
Netzer 2009). Since the dependence on distance to the light emit- 
ter is an inverse square law, like for gravity, the effect of radiation 
pressure support is to reduce the apparently measured mass of the 
black hole, based on BLR kinematics. Recently, it has been un- 
der de bate, how strong this effect was dMarconi et al.ll2008ll2009l ; 
iNetze r 2009; Netzer & Marziani 2010|), especially in the case of 
NLSl galaxies. It is interesting to ask in this context, which cloud 
orbits would be stable against small perturbations, and also, which 
ones are compatible with the spectropolarimetric results. In the fol- 
lowing we perform such an analysis. 
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Figure 1. Dynamical equilibrium column density over luminosity in Ed- 
dington units against rotation velocity in Kepler units for an isotropic light 
source. The red, dashed part of the line corresponds to a maximum of the 
effective potential for reasonable choices of the parameter i (.s > 1) that 
characterises the pressure profile. Stable orbits are still found in this case, 
but are highly eccentric, and are found above the red line. The yellow, dot- 
ted part corresponds to a stable minimum for certain values of s. The solid, 
black region is always a minimum of the effective potential, provided i < 3, 
and therefore allows for orbits with low eccentricity, which scatter around 
the line. See text for more details. 



We present the cloud orbit analysis in sect.|2l discuss our find- 
ing in sect.[3]and summarise our conclusions in sect.|4] 



2 DYNAMICAL STABILITY ANALYSIS 

We carry out the analysis for both, an isotropic light source, and a 
light source with a cos Q dependence for the luminosity, where is 
the polar angle. The latter case should be more realistic, as it is the 
expected angular distribution for the luminosity of a thin accretion 
disk. 



2.1 Isotropic light source 

For optically thick, spherical clouds of constant mass m and central 
column density A', without internal structure, the force equation in 
spherical polar coordinates reads: 



GMbhot / 3/ 



R^ yiajN 



(1) 



where C7x is the Thomson cross section, / the luminosity in Edding- 
ton units, V the rotational velocity in Kepler units and Mbh the 
mass of the black hole. 

Dynamical equilibrium, corresponding to circular orbits, is 
reached for a column density of: 

3Z 



N-- 



2(7t(1-V2) 



(2) 



For a large range of rotational velocities V , this is of order 
10^'' / cm^^ (compare Fig.Q}. 

To assess dynamical sta bility, we now consider perturbations 
to circular orbits. Following iNetzen d2008l) . we assume a confin- 
ing inter-cloud medium with a pressure profile of p{R) oc R^'' . 
The sound crossing time through a BLR cloud is of order weeks, 
whereas the orbiting period is many years. We therefore assume 



the clouds to adjust their radius instantaneously to any change in 
external pressure. A change in the clouds' cross section alters the 
amount of radiation received and hence the radiation pressure sup- 
port. Clouds will be dynamically stable, if an increase (decrease) 
in orbital radius R leads to a net inward (outward) force, i.e. a po- 
sition above (below) the equilibrium line in Fig. [T] We proceed to 
calculate the perturbed locations in the N/l — K-diagram. 

We first consider a spherically symmetric cloud initially in a 
circular orbit at distance R to the black hole. Since the temperature 
is kept constant by photoionisation, a change in distance by A/f > 
will result in a change of the particle density by n{R + AR) = 
n{R)[{R + AR)/R]^'\ This results in a change of the column den- 
sity: 

N{R + AR)=N{R){1+AR/Ry^'^^. (3) 

The corresponding change in orbital velocity for angular momen- 
tum conserving perturbations is: 

-1/2 



V {R + AR) = V {R){1 + AR/R) 

The perturbed cloud will receive a restoring force, if: 



N(R + AR) > 
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(5) 



It may be shown by a few lines of algebra that this condition may 
be fulfilled only if 
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(6) 



The analogous derivation for tJi < leads to the same result. If 
the rotational velocity is below this value, positive perturbations to 
R lead to ejection. Negative ones lead to highly eccentric orbits. 
We show this by considering the effective potential Veff- Inserting 
eqs. ([3} and Q into eq. (TJ, and defining Rq to be a fiducial distance 
to the black hole where eq.|2]holds, with N{R()) = No, V{R()) = Vq 
and X = R/R(), we derive: 
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s = 3/2 
1 5^3/2 

The effective potential is displayed in Fig. |2] In agreement with 
the preceding discussion, there is an extremum at x = 1, whose 
character depends on Vq: For Vq > K: (compare eq.|6), it is a mini- 
mum. Stable bound orbits with small radial motions may be found 
in this case. For Vq < Vc, the extremum at jc = 1 is a maximum. 
In this case, their exists a minimum further in. The equations can 
be solved easily analytically for i = 0, 3/2 and 3. For .s = 3/2, the 
second extremum is at: 
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This extremum is a minimum for Vq < Vc . We have verified numer- 
ically that for all 1 < .s < 3, Xe2 is very close to the value of the 
s = 3/2-case, for Vq < Vc. Again, stable bound solutions may be 
found. If a cloud is very deep in that potential well, the orbits are 
close to circular, and dominated by rotation. The average column 
density and rotation velocity for such an orbit are much higher than 
the values at Rq. Such orbits are therefore effectively on the stable, 
solid black part of the equilibrium curve in Fig. [T] Orbits with to- 
tal energy close to the potential energy at the maximum will follow 
highly eccentric orbits, with predominantly radial kinematics. For a 



Stability of Orbits in the Broad Line Region 3 
Potential - isotropic, Vo=0.1 








^ ' 1 







- - ^ 


-_ 




/ <■■' 


■■"v-.... 


10 


' ' /'" 


\ ''■.. 


-- 






\ 




20 


! f 


\ 


'■--..^ i 




■\ / f 






30 


11 A 


\ 


J 




y .'" 




; 


40 


r • / 




i 




'. 1 




t 


50 


_ \f 




i 



-2-10 1 2 

log (X) 

Potential - isotropic, Vo=0.9 






0.0 



-0.5 



-1.0 - 
-1.0 



-0.5 



0.0 

log (X) 



0.5 



1.0 



Figure 2. Normalised effective potential against normalised distance be- 
tween cloud and black hole for a rotation velocity of Vq = 0.1 (top) and 
Vq = 0.9 (bottom) normalised to the Kepler value at R = i?o (.v = 1) for 
the case of isotropic illumination. In each case, the solid green (dotted red, 
dashed yellow) line is for s= \ (2,3). 



significant fraction of their orbital period, they have indeed low col- 
umn densities and rotational velocities, corresponding to the region 
above the red dashed part of the line in Fig.[T] 



2.2 Anisotropic light source 

The luminosity in Eddington units for a geometrically thin accre- 
tion disk is a function of the polar angle and is given by 2/| cos 6 1, 
where / is the total luminosity of the source in Eddington units. 
The force is still central and therefore, angular momentum is con- 
served. The orbits are planar. Any particular orbital plane may be 
characterised by a polar angle do- We define the direction of the 
maximum elevation above the equatorial plane to have an azimuth 
of (^ = 0. A cloud will then have maximum radiation pressure sup- 
port at (^ = and (j) = 7t, and no radiation pressure support at 
(p = ;r/2 and (j> = 3;r/2, when passing the equatorial plane. Ge- 
ometrical considerations result in: cos 6 = cos (j) cos 6o- Using this, 
and eqs (TJ, ([3} and ([4}, results in the following total force equation: 

^_ GMBHin f 31\cos(j)coseo\ ; , \ 
^^^[ CTTiVox-2V3 +^«" "V' ^^ 

where x = R/Rq, R{(I) = k/2) = Rq, N((t) = Tl/2) = Nq. The az- 
imuthal velocity at (^ = Tt/2 in Kepler units is denoted by Vq. Be- 
cause of the (^-dependence, the force is no longer conservative. We 
therefore expect the total energy and hence also the major axis of 
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Anisotropic light source: orbit parameter 




Figure 3. Force-free lines for the case of the anisotropic light source: nor- 
malised distance between cloud and black hole versus a\ cos (^ \ for different 
values of the pressure distribution index .v. The force is outwai'd to the right 
of each curve, and inward on the left side. Around the lower branch of each 
line, the force is always restoring. Stable orbits are therefore expected in this 
region of the parameter space. If the cloud has low enough column density, 
a too large part of the orbit will be right of the force-free line, where the 
force is always outward, leading to ejection of the cloud. 



the orbits to grow or diminish, depending on whether the cloud 
moves predominantly with or against the radiation flux. 

First, we consider the force free locations, which are useful to 
understand the general structure of the allowed orbits. In the stan- 
dard treatment of the Kepler problem, bound orbits require the ex- 
istence of a local minimum of the effective potential. The elliptic 
orbits are then oscillating around this minimum. The forces in this 
case vanish, of course, at the minimum. Away from the minimum, 
there is an effective restoring force. If the total energy is too large, 
the object can reach a region where the effective force is away from 
the equilibrium point, and may escape. A similar analysis may be 
done for the present problem: Making the definition: 

3/COS0O 



(10) 

we find the condition for force-free lines: 

a|cos0|=x2^/-^(l-Vox-'). (11) 

A force-free line in the orbit plane is defined by specifying the pa- 
rameters a and s. Since there is no radiation force in the equato- 
rial plane, force-free lines require Vq = 1 • We show x as a function 
of a I cos (^ I for different values of s in Fig. [5] Each curve is char- 
acterised by a rightmost value, a|cos(|)| = Omax. where the upper 
branches with negative slopes and the lower branches with posi- 
tive slopes meet. The force is restoring in the vicinity of the lower 
branches with positive slope. Therefore, if a cos </i < OmaxV^ , we 
expect bound orbits to exist. We do not necessarily expect bound 
orbits if a cos exceeds this value for some fraction of the orbit. 
From eq. l IlU . we find: 



2s \ 1 



1 



=^ + 1 



(12) 



For 5 = 2 (1,3), flmax evaluates to 0.20 (0.33, 0.15). Consequently, 
there is a 0o and /-dependent critical column density, above which 
bound orbits should be found: 

COS0n 0.33 



A'o>7x 10^^ cm" 



0.1 0.5 a„ 



(13) 



Similarly as in the isotropic case, we also find bound, but 
highly eccentric orbits for low column densities and rotation veloci- 
ties: We show this by numerical integration of some example orbits 
(Fig|4l(. Our reference in the following is a Cartesian coordinate 
system with coordinates X and Y defined in the orbital plane. The 
K-axis is taken to be the intersection between the orbital and the 
equatorial plane. We use s = 2 and start the clouds at = k/2, i.e. 
where the orbital plane meets the equatorial plane (positive K-axis 
in Fig|4ll, with certain values for a and Vq. For (a, Vq) = (0-1, 1-1), 
we expect a stable bound orbit despite the super-Keplerian initial 
velocity, since a is in the stable regime. This is indeed the case 
(Fig. 151 top left). For (a,Vb) = (1,1), we expect ejection, since a 
is in the unstable regime and the comparatively high velocity cor- 
responds to a positive radial perturbation. Again, this is what we 
find (Fig.|4] top middle). For (a,Vb) = (1,0.2), we expect highly 
eccentric orbits. This is confirmed by the numerical integration 
(Fig- ID top right). Here, the non-conservative nature of the poten- 
tial is most apparent: The cloud gains energy, when moving out- 
wards. But since it advances in azimuth, it does not get back the 
same amount on the way inwards. However, on average over many 
orbits, the contributions cancel each other. 

In order to obtain kinematic information, we also recorded the 
emission measure, e, of the cloud into a certain direction during the 
integration, together with the current position and velocity. At each 
time /, it is calculated as: 

e{t)=A{(l>)AtR-''l^-^, 

where A((^) is the fraction of the illuminated cloud surface seen by 
the observer. We take only the side facing the AGN as line emitting 
region, as appropriate for an optically thick line. The current time 
step interval is denoted by At, and the radial dependence is due to 
the change of the cloud size and radiation flux with distance to the 
centre. We place two observers at large negative X and Y values, 
respectively. At each timestep, we project the clouds velocity along 
the respective lines-of-sight. e{t) At is then added to the correspond- 
ing bins of line-of-sight velocity and transverse position in order to 
create a two-dimensional emission-weighted histogram. 

For the bound orbits, these velocity resolved emission- 
weighted histograms are shown in Fig|4] A priori, one might expect 
a broader signal for the low a (high column density) case with well 
separated emission peaks near the positive and negative Kepler ve- 
locity, and a narrower signal for the high a (low column density) 
case. This is indeed what we find. For (a, Kq) = (0.1, 1.1), the peak 
of the emission is at ± f« 50 per cent of the Keplerian velocity 
at R(^. For (a, Vg) = (1,0.2), the orbits get very anisotropic in real 
space as well as in velocity space. An observer in the — X-direction 
would see two peaks at around ±0.15 times the Kepler velocity 
at i?o, close to the initially imposed one. Here, the outermost lo- 
cations of the orbit dominate the emission due to the longer time 
spent there. An observer in the — F-direction would see the two 
peaks at around 1 .2 times the Kepler value, because from this point 
of view, the orbits are much nan^ower. If there was an ensemble of 
such clouds with the angular momentum vector randomly rotated 
around the symmetry axis of the system, but otherwise identical, 
the signal at low velocity would dominate, as the emissivity at the 
slower peak is about three times higher. For an observer who would 
see the orbital plane at some inclination, the apparent velocities 
would be somewhat below these values. Since spectropolarimetric 
observations are able to resolve the BLR in many objects, we have 
also separated the emission that comes from the positive part of 
the transverse axis from the one of the opposite side. As one might 
have expected, the emission from the two sides is well separated in 
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X-observer: velocity histogram 




Figure 4. Position and kinematic data for three example orbits for the anisotropic case. Radiation pressure support is maximal on the X-axis, and zero on the 
Y-axis. We have integrated for 100 orbits, for the bound cases. The simulations are converged with regard to the timestep. Top row: orbit shapes in orbit plane. 
The orbits are similar to precessing Kepler-ellipses. Because of the non-conservative nature of the potential, they grow and shrink periodically, with a period of 
Tl. This leads to Rosetta shapes, which are compressed in the direction of maximum radiation support. The set of four plots on the lower left shows kinematic 
information for the (a,Vo) = (0.1, 1.1) simulation. Here we assume that a cloud emits radiation on the illuminated side facing the centre. We calculate the 
fraction of that surface seen by respective observers. The emission is further proportional to the current cloud surface and the time interval. The top row 
relates to an observer in — X-direction (X-observer), the bottom one to an observer in — F-direction (Y-observer). The left column shows position-velocity 
diagrams, the right one the emission at a given velocity (black) and separately for negative (red) and positive (blue) values of the respective complementary 
coordinate. The emission of the two peaks at about half the Kepler velocity is easily discerned. The lower right set of four plots shows the same for the run 
with (a, Vb) = (1,0.2). The emission occurs predominantly at smaller velocities. In velocity space, emission from the sides with positive coordinates still has 
little overlap with the emission from the other sides. Colour scales are cut at reasonable values, as indicated at the individual colour bars. Velocity is given in 
units of the local Kepler velocity at the initial cloud position Rq. 



velocity space, in both cases. Remarkably, for the emission from a 
given side, the (a, Kq) = (1,0.2) cloud shows about ten percent of 
the peak emission of the opposite side at the location of the peak of 
the opposite side. In contrast, for the cloud with (a, Vq ) = (0. 1 , 1 . 1 ) , 
the emission of a given side drops to zero at the peak of the other 
side. 



3 DISCUSSION 

We have shown that pressure confined clouds at any sub-Keplerian 
rotational velocity may exist in stable dynamical equilibrium in the 
BLR. Given the significant evidence for a gravitationally bound, 
flattened, but still of considerable thickness, and disk-like geometry 
for the BLR (compare sect. [TJ, it is reasonable to require a stable 
dynamical equilibrium for the line-emitting clouds. 

An essential ingredient for the model is the behaviour of 
the pressure of the inter-cloud material with radius: A reason- 



able assumption for the inter-cloud comp onent is an Advec- 
tion Domin ated Accretion Flow (ADAF, e.g. iNaravan & Yill 19941 : 
lYuan. Ma, & Naravan 2008). A nearly hydrostatic solution is in- 
cluded in the ADAF models for the limit of low accretion rates. 
For this type of solutions, the power law index for the pressure s 
is between two and thr ee. Outflow solutions ha ve also been con- 
sidered in the literature. [Konigl & Kartid ( 1 19941) find 1 < .y < 1.5. 
Therefore, . ; should be between one and t hree (similar results are 
obtained bv lRees, Netzer. & Ferlandl 19891) . 

We have considered an isotropic and an anisotropic light 
source as appropriate for an accretion disk. The former is usually 
implied in the literature. For the isotropic case, the force is cen- 
tral (i.e. conserves angular momentum) and conservative. We find 
an equilibrium relation between column density, luminosity of the 
AGN and rotational velocity. We show by direct stability analysis 
that only for the part of that relation with high rotational velocities 
a cloud would encounter a restoring force for small radial perturba- 
tions. However, by analysis of the effective potential, we show that 
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stable orbits may be found also for low rotational velocities and col- 
umn densities slightly above the equilibrium curve given by eq. ^. 
However, the character of these orbits changes: while high column 
densities (corresponding to close to Keplerian rotation) allow for 
even circular orbits, low column density clouds require highly ec- 
centric ones. 

For anisotropic illumination, the force is still central, and 
therefore angular momentum conserving, but it is no longer conser- 
vative. Therefore, as the orbits precess, there is now the additional 
feature that they gain and loose total energy, which is exchanged 
with the radiation field. This is evident from the periodic change of 
the major axis of the orbits. We show that the time averaged signal 
of such a cloud would also be very anisotropic. If one would con- 
sider many such clouds with randomly rotated angular momentum 
vector around the axis of symmetry, the stronger emission from fur- 
ther out would dominate. In our example, this produces a peak in 
emission at a small fraction of the Keplerian velocity with a similar 
FWHM. More important, the emission contributions from the two 
sides of the accretion disk are kinematically clearly distinct. This 
would therefore probably still be compatible with the spectropo- 
larimetric results. More detailed comparisons using realistic cloud 
samples have to be done in order to decide if the dispersion would 
be low enough to fit with all the available data. This is however 
beyond the scope of this article. 

Cloud column density and rotation velocity in the anisotropic 
case are still related by eq. ^ (also Fig. [T), when orbit averaged 
values are used. Therefore, it is well possible to observe arbi- 
trarily low rotation velocities for cloud col umn densities of orde r 
10^-'(//0.1) cm"-, confirming the results of iMarconi et al.l l llOOa) . 

Interestingly, it has turned out that for the anisotropic case, 
the orbits are less extended in the direction of the strongest radia- 
tive force. This might first appear counter-intuitive, but is readily 
explained if one considers that the radiative force adds energy to 
the orbit, as long as the motion is outward, but brakes down the 
cloud, when it moves inward. If the major orbit axis coincides with 
the direction of maximum radiation force, the contributions nearly 
cancel, whereas one may get positive changes when the major axis 
has advanced past that direction. Energy losses are expected, when 
the major axis has not yet reached the maximum force direction. 

If the column densities are not too high, the radiation field 
favours a disk configuration: Consider an initially isotropic ensem- 
ble of clouds with relatively low column densities, comparable to 
the one given in eq. ( 113b . and a broad distribution of angular mo- 
menta. Due to angular momentum conservation, each cloud orbits 
the black hole in a particular orbital plane on an elliptical orbit. 
Once an anisotropic AGN is switched on, we may assign an a 
value to each cloud. Because a <x cos(0o), clouds at low polar angle 
(closer to the symmetry axis of the emission of the central accretion 
disk) have greater a values. If this a value is too large, the high an- 
gular momentum clouds will be ejected. For the remaining clouds 
at low polar angle, the major axis of their orbits will shrink when 
it points towards higher latitudes on their course of precession. In 
Fig|4]we have demonstrated this effect. The result will be a disk- 
like BLR. For column densities much higher than the one given in 
eq. ( I13t , the orbits are less affected by radiation pressure. For such 
clouds, the BLR would therefore not be constrained to have a disk 
or other geometry by these dynamic considerations. It may or may 
not be in a disk configuration for other reasons. The critical col- 
umn density is of the order 10""* cm"- for common luminosities 
of ten per cent of the Eddington value, which is rather large com- 
pared to observational constraints (compare references in sect. [TJ. 



One might therefore generally expect this effect to be significant in 
many BLRs. 

The detailed mixture of orbits is of course very hard to pre- 
dict. Since the cloud mass is unimportant for the acceleration of 
the cloud, we expect clouds of a wide range of masses, with 
the column density adjusting according to the cloud's position in 
phase space. Some B LR clouds might have been bom in situ (e.g. 
IPerrv & Dvsonll 19851) . If the clouds would have come from further 
out, an interaction would be required to reach the bound orbits. This 
might favour eccentric orbits. The mixture of orbits will determine 
the observed velocity structure. 

Pressure confined spherical clouds as used in this analysis, 
suffer from shearing by differential radiative forces due t o the vary- 
ing c olumn density from the cloud's rim to it's centre l IMathewsl 
119861) . For our parameters one would expect complete disruption 
after about a hundredth of an orbital period. This issue is com mon 
to this class of cloud model (compare e.g. iRees et al.lll989l) . and 
has not been solved so far. Possible ideas to stabilise the cl ouds in- 
clude magnetic fields and a more favourable geometry (e.g. lNetzeg 
l2008l) . 



4 CONCLUSIONS 

We have shown that pressure confined clouds may rotate stably 
on bound orbits near the dynamical equilibrium between radia- 
tion, centrifugal and gravitational forces at all sub-Keplerian ro- 
tational velocities. This is true for isotropic illumination, as well 
as for the case where the radiation flux is correlated with the po- 
lar angle. While angular momentum is conserved in both cases, 
energy is not conserved for anisotropic illumination. This leads to 
Rosetta orbits that extend less in the direction of maximum radi- 
ation force. An intrinsically isotropic low column density cloud 
system would therefore become less extended in the polar direc- 
tions, when an anisotropic AGN would be switched on, and conse- 
quently appear disk-like. We show that it is possible to find clouds 
of low and high rotational velocity with well separated peaks in 
the spatially resolved emission spectra as a function of velocity, 
as required by spectropolarimetric BLR data. These findings con- 
firm the idea that significant corrections of black hole masses due 
to radiative forces ar e possible in certain objects, as proposed by 
JMarconi et alj2008t) . 
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